Josephson currents through spin-active interfaces 
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The Josephson coupling of two isotropic s-wave superconductors through a small, magnetically 
active junction is studied. This is done as a function of junction transparency and of the degree of 
spin-mixing occurring in the barrier. In the tunneling limit, the critical current shows an anoma- 
lous T~ l temperature dependence at low temperatures and for certain magnetic realizations of the 
junction. The behavior of the Josephson current is governed by Andreev bound states appearing 
within the superconducting gap, A, and the position of these states in energy is tunable with the 
magnetic properties of the barrier. This study is done using the equilibrium part of the quasiclassi- 
cal Zaitsev-Millis-Rainer-Sauls boundary condition for spin-active interfaces and a general solution 
of the boundary condition is found. This solution is a generalization of the one recently presented 
by Eschrig [M. Eschrig, Phys. Rev B 61, 9061 (2000)] for spin-conserving interfaces and allows 
an effective treatment of the problem of a superconductor in proximity to a magnetically active 
material. 

PACS numbers: 74.50.+r, 74.25.Ha, 74.80.-g, 74.80.Fp 



I. INTRODUCTION 

If a supercaaductor is exposed to magnetically ac- 
tive impuritiead'EI or materialsQ the superconducting state 
is modified. Josephson coupling two superconductors 
through a magnetically actiwe barrier may lead to what 
is known as a "7r"-junctionB, a junction for which the 
ground state has an internal phase shift of 7r between the 
superconductors across the barrier. If the barrier is ex- 
tended to an S/F/S-structure, a ferromagnetic (F) layer 
sandwiched between two superconductors (S), the criti- 
cal current will oscillate as the thickness of F is variedo. 
Additionally, the critical current will also depend on the 
strength of the exchange field in F. The principal rea- 
son for this strong dependence of junction properties is a 
drastic modification of the local superconducting density 
of states in the contact region to the F-layerD. Scat- 
tering of a magnetically active surface or transmission 
through a ferromagnetic barrier leads to a depairing of 
the Cooper-pairs and the creation of surface or layer An- 
dreev states at energies within the superconducting gap. 
There have been an extensive experimental effort to ex- 
plore the physics |-aJpove by tunneling through magnetic 
insulator barrigrsEra, by probing the proximity effect in 
S/F-structuresQ and constructing S/F-multilayersB (see 
also references therein). The problem at hand is quite 
formidable since the strength of the exchange energy is 
for most ferromagnetic materials, like Ni, Co and Fe, a 
sizable part of the Fermi energy (~ eV) while supercon- 
ductivity lives on a much smaller energy scale (~ meV). 
To overcome the difference in energy scales the ferromag- 
netic layer must be extremely small (< nm) and only re- 
cently have supercurrents been reported in S/F/S junc- 
tions by Veretennikov et a£l using weak ferromagnetic 
alloys for the F-layer. 

To efficiently model the S/F/S-junction there are two 
main routes of approach. The first is to assume an ex- 



tension of a ferromagnetic metal, now characterized by 
a length and an exchange field, separating the two su- 
perconductors. jiiZlthin this approach both critical cur- 
rent oscillationstrt2l and the effecjUpf the exchange field 
on the Andreev bound statesE3li3 have been studied. 
The limitation of the approach is that it is restricted 
to small exchange fields, i.e. fields that are compa- 
rable to the superconducting gap. An alternative ap- 
proach is to treat the ferromagnetic part as a partially 
transparent barriac. [Which transmits the two spin pro- 
jections differentlya'Q'Ejta. Using Bogoliubov-deGennes 
equations and a WKB-approach for the ferromagnetic 
barrierO, JosephspHpGurrent-phase relational^ and quasi- 
particle tunneling! 1 3Q have been studied for both con- 
ventional s-wave and unconventional d-wave supercon- 
ductors. 

In this Paper, I follow the second path making use 
of the quasiclassical theory appropriate for describing 
low-energy phenomena like superconductivity for which 
considered energies are small compared to the Fermi 
energy Ep [ |lS|-|2C|]. Ferromagnetism will enter as a 
boundary problem for the quasiclassical Green's func- 
tions g{pf, R s 'i e) at a semi-transparent interface separat- 
ing two conventional s-wave superconductors. There is no 
general restriction of validity of the present work to con- 
ventional superconductivity. The physics revealed in the 
simplest system proves to be quite rich without addjiia 
properties related to an unconventional pairing stateEjtLZl 
and therefore s-wave superconductivity in proximity to 
ferromagnetism should be studied in its own right. In sec- 
tion [n|, a general solution of the quasiclassicaLboundary 
condition, as posed by Millis, Rainer and Saulalj, is given 
for equilibrium Green's functions. In section III the local 



density of states in proximity to a ferromagnetic insulator 
is discussed. This is an important step in understanding 
the Josephson coupling between two s-wave superconduc- 
tors studied as function of a simple phcnomcnological two 
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spin-band scattering S'-matrix. Section IV is devoted to 
the study of the Josephson coupling and maps out regions 
where the junction is in a normal "0" -state and where it 
switches to the "7r" -state. 



II. BOUNDARY CONDITIONS FOR 
QUASICLASSICAL PROJECTORS AT 
SPIN-ACTIVE INTERFACES 

Surfaces and interfaces involve energies of order Ep 
which in quasiclassical theorytS are integrated out at 
the onset. This means that boundary conditions for the 
quasiclassical Green's function at surfaces and interfaces 
must be posed for the full Green's function satisfying 
the Gor'kov equation. Resulting boundary conditions 
hava-then to be energy integrated into their quasiclassical 
formEy. Physical properties of an interface can then be 
accounted for by a suitably chosen scattering S-matrix. 



A boundary condition for partially transmitting inter- 
faces was first rdpived by Zaitsevt3 and, independently, 
by KieselmannE2l. The boundary condition was later gen- 
eralized by Millis, Rainer and Sauls (MRS) to include 
spin-active interfaces!!^!, i.e. interfaces which transmit 
and reflect quasiparticles differently depending on the 
spin projection. _ 

Recently, EschrigcJ used a projector method to solve 
Zaitsev's boundary condition in general. The projectors, 
Va^ , introduced relate to the quasiclassical Green's func- 
tion as = —in(Pjp — V^). The superscripts (I) 
label the side of the interface, the subscripts ± are direc- 
tional indices and finally the " hacek" denote the Keldysh- 
matrix structure of the Green's function. For a full ac- 
count on the quasiclassical projectors Va^ and their pa- 
rameterization, I refer the reader to Eschrig's original 
paperO and in Appendix ^ I give a brief review of ele- 
ments of quasiclassical theory used in this paper. Written 
in projectors Va ^ equations (63-66) of MRS reads 
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Here, the non-commutative ^-product is a usual matrix 
product and a folding of internal energies. The solution 
of the system of equations ([!]) is facilitated by a conve- 
nient parameterization of V a in terms of four coherence 
functions -f R ,j R . 



x K ,x K . 



"f A , j A and two distribution functions 
These six functions are 2x2 spin matrices and 



superscripts (R,A,K) stand for Retarded, Advanced and 
Keldysh. The set of functions above obey Riccati-likc 
equations that are easier to .handle than the original 
quasiclassical matrix equations. Especially, they fulfill 
certain stability criteria when integrated for along tra- 
jectories x. The functions 7' R ,7" 4 and x K are bounded 
when integrating the along a trajectory (vf(p) ■ x > 0) 
and functions j R , ^ A and x K are bounded integrating in 



the opposite direction (f/(p) ■ x < 0) [ 25|. Restating 



this in context of the interface problem we can always 
integrate up to the barrier obtaining 7^,7^, x K along 



trajectories with Vf 
p 1 and p 2 in Figure 



t)-x ; 
Simi 



> 0, i.e. along the trajectories 



Similarly, 7 , 7 and x are inte- 



grated stably along Vf(p) ■ x < 0, i.e. p 1 and p 2 in Figure 
lL For constructing the Greens's function at the surface 



that fulfills eqs. ([!]) we still need the "scattered" func- 
tions T R , T A and X K with v f (p) ■ x < and T R , T A and 
X K with Vf(p) x > 0. These functions have to be solved 
for from eqs. ([!]). Using the functions {j R ,J R , T R ,T R ), 
the Retarded projectors are constructed as 
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FIG. 1. A schematic picture of the in-scattering trajecto- 
ries (p 1 ,p 2 ) connected over an interface barrier parameterized 
by an S-matrix to the out-scattering ones (p ,p ) 
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and after substitution into eqs. (P one finds after some 
straightforward algebra that the scattered-out functions 
can be expressed solely by scattering-in functions as 



rf 



(2) 



The generalized reflection coefficients TZ R are defined as 



lZf r — Snp 21 1 « [Snp 21 1 - S12P22 1 ] 
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and corresponding transmission coefficients T R p = 1 
TZ§ 1 , The functions p R and are given as p R 



(3) 



5y - 7 f ®<%7.f and p£ = 5?. - 7 f «5t.7f . The scat- 
tered coherence functions on side 2 are given by inter- 
changing the side index 1 «-> 2. Advanced functions 
are related to the retarded ones by general symmetry 
g A = T3(<7 fi )^T3. The similarity in form of the final re- 
sult given in equation (||) to the solution of a scattering 
problem is not a coincidence as the boundary condition 
for the quasiclassical Green's ftuaction can also be solved 
by a direct scattering approached. 

So far no reference to the form of the scattering 5- 
matrix has been made. 5 is a scalar in Kcldysh space 
and a matrix S in particle-hole space, spanned by Pauli- 
matrices fj, with the form S = 5(1 + r 3 )/2 + 5(1 — t 3 )/2 
where S(p..) = 5 tr (— p,,) [ ^3). For spin-active inter- 
faces the different components of the 5-matrix, 5y in 
(||) above, are 2x2 spin matrices. To proceed further a 
specific 5-matrix is chosen to model the magnetic barrier 
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where aj notes the Pauli-matrices spanning spin space 
and parameters (t, r) are the usual transmission and re- 
flection coefficients. The S-matrix (^) is one of the sim- 
plest choices that allows a variable degree of spin mixing 
at the interface and the spin mixing is parameterized by 
the spin-mixing angle O. By this construction 5 only vio- 
lates spin conservation, i.e. it does not commute with the 
quasiparticle spin operator <x. The angle O will be con- 
sidered as a phcnomenological parameter independent of 
the trajectory direction in this paper, but as shown by 
Tokuyasu et al in the appendix of Ref. [ 27j one can 
relate to the microscopic properties of the magnetic 
barrier. In particular, in Ref. [ 2jJ an S-matrix is con- 
structed for a magnetically ordered insulating barrier and 
it is found, as expected, that 5 depends on the quasipar- 
ticle momentum projection parallel to the interface and 
on material parameters describing the barrier such as the 
average band gap, E g , the internal exchange field, hi, and 
its orientation (1. 

Only the simplest case of an isotropic s-wave supercon- 
ductor will be considered in this paper using weak cou- 
pling BCS theory. For this, assuming a constant order 
parameter in space A(x) = A, the retarded coherence 



functions are j R — 70*02 and 7^ = 10270 with 70 = : 
-A/{e R + in), 70 = A*/(e R + in), Q = ^|A| 2 - {e R f 
and e R = e + i8. If, on the other hand, the effect of 
proximity to a magnetic material oa the superconductor 
is of interest the Riccati equationscil 



i v r V R ~{ R + 2e R -/ R = 7 R A R j R + ~/ R - 7 
i v r V R j R ~ 2e R j R = ^y R A R j R + E R j R - ~ R r R 



R T, R - A R 
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have to be solved together with a self consistent determi- 
nation of the order parameter A(x) and of the impurity 
self energy S(x;e) as described in the appendix. In eq. 
(||), A R = A + T, R d and A R = A + t R d are the impu- 
rity renormalized gaps, while H R and T, R are diagonal in 
Nambu space, and include both the impurity self ener- 
gies Ej and the mean fields S m . All functions entering 
are spin matrices, and for the problems that will be con- 
sidered in this paper it is sufficient to parameterize the 
matrices by two components as e.g. ( 7o + 7303)102 and 
i<72(7o — 7303)- The bulk values for j R and 7 , given 
above, serve in this case as initial values when integrat- 
ing eqs. (§). 



III. ANDREEV BOUND STATES AT AN 
IMPENETRABLE MAGNETIC BARRIER. 

As a first poiBt, we return to the half-space model 
of Tokuyasu et a£l, a semi-infinite BCS-superconductor 
bounded by a magnetic insulator. The scattering off the 
insulator is assumed to be specular with a phase shift ac- 
quired differently at reflection for spin-up and spin-down 
quasiparticles. Using the 5-matrix above, the coherence 
functions scattered off the magnetic insulator are given 
directly by the incoming ones as T = cxp (2003)70202 and 
r = i(J2 7 o exp (i0<73). Using the information of the scat- 
tered coherence functions the Green's function is given 
for a trajectory with n • py > as 



-in N 1 



(1 + j R f R ) 2 7 « 
-2f R -(l + f R j R ) 



where 



N R = 



(1 - 7 ii f fi )- 1 

(l-fV*) -1 



(6) 



(7) 



If the trajectory is n ■ < 0, the Green's function is 
simply given by interchanging 7 — > T and T — > 7. The 
effect of the spin mixing is perhaps best seen in the spin 
and angle-resolved density of states (DOS) right at the 
barrier. On a trajectory with fi ■ py > this quantity is 
given by the imaginary part of g R {pf, s), th e upper left 
component of (||) as described in eq. ( A10Q . Assuming 
a constant order parameter up to the interface, the local 
DOS at the interface for spin-up quasiparticles, Aj, can 
be written as 
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FIG. 2. Tunneling density of states for spin-up quasiparticles at an impenetrable magnetically active insulator. The 
superconductor is a conventional BCS superconductor with a self consistently determined order parameter. In the main 
panels is varied from to n in steps of n/4 going top to bottom and the curves are shifted for clarity. In the left 
figure the impurity scattering is in the Born limit while in the right panel the scattering is taken in the unitary limit. In 
both cases the scattering rate is chosen to 0.01 in units of 2nT c . In the left inset in the left figure the DOS is plotted for 
Q — ty but with varying impurity scattering. The dependence of the zero-energy peak width with y/T is plotted in the 
right inset of the left panel. In the panel to the right the impurity scattering is in the unitary limit. As expected the 
broadening of the Andreev states is very much suppressed. 
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For spin-down quasiparticles N± reads the same after sub- 
stitution 0^—0. This density of states has Andreev 
bound states inside the gap, i.e. for |e| < A. These 
states are located at £b,\(i) = ±Acos(0/2), with +(— ) 
for the spin-up (-down) branch. It is notable that the 
DOS (||) has exactly the same form as the DOS calcu- 
lated for a S/F/S weak link, for which the angle is 
shown to depend on the thickness of and exchange field 
in the F-layerEa. 

The existence of bound states will lead to a reduc- 
tion of the order parameter amplitude in the vicinity of 
the magnetic interface. This is seen in Figure 2 in Ref. 
[ Enj. The pairbreaking occurs gradually as the bound 
state on either spin branch is tuned towards £& = with 
— > 7r. At = 7T the order parameter is totally re- 
duced at the interface and recovers to its bulk value over 
a distance of order £ G , the zero temperature coherence 
length in the superconductor. It turns out that the po- 
sition in energy of the surface states, is not to sensi- 
tive to a spatial variation in the order parameter. The 
Andreev bound states would be delta peaks if the super- 
conductor had an infinite mean free path. The presence 
of impurities in the bulk give rise to a finite life time r 
which broadens the Andreev peak. The broadening of 
the Andreev states is sensitive to the scattering strength 
of the impurities. In the main panels of the figures || 
the DOS is shown for different values of spin- mixing 0. 
The surface order parameter is self consistently deter- 
mined for each spin-mixing angle. The scattering rate 
is T = 1/2t = 0.01 in units of 2irT c and corresponds to 



a mean free path £ m f p = 50 £ D . From the literature on 
d-wave superconductors it is known that surface bound 
states are broadened by impurity scattering as ~ \J Ar 
in the Born limito. In the two insets in figure || the 
dependence of the width of the zero energy peak with 
scattering rate for = 7r is shown. For small scatter- 
ing rates indeed the V AT-dependence is recovered. If, 
on the other hand, the scattering strength of the impu- 
rities is in the strong scattering limit the broadening of 
the Andreev bound states will be exponentially small, 
~ \J Ar exp[— A/r] [ ^9|. The effect of unitary scatterers 
is shown in the figure to the right in Fig. 0. 



IV. JOSEPHSON CURRENT-PHASE RELATION 
AND THE ENERGY STATE OF THE JUNCTION. 

Next, let us consider the Josephson current-phase rela- 
tion through a magnetically active point contact. Assum- 
ing a point contact allows several simplifications. Effects 
of the contact itself on superconductivity i.e. the order 
parameter profile may be disregarded. This holds true if 
the contact radius is taken much smaller than the super- 
conducting coherence length. Furthermore, spin-neutral 
surface scattering alone does not affect an isotropic s- 
wave superconductor and thus bulk values of the coher- 
ence functions 7, 7 can be used for the in-scattering ones 
in the boundary condition (|2|). An additional advantage 
of the point contact condition is that the results will not 
depend on the presence of non-magnetic bulk impurities 
as the current through the point contact depends only on 
bulk coherence functions. On the other hand, the point 
contact itself is fully described by its transmission t and 
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degree of spin mixing, 0. The Josephson current through 
the contact is calculated as a function of the phase dif- 
ference, <f>, between two superconductors by the current 
formula 



j{4>) = eN f j —Tr(v f T 3 gf(p f ,0; e))p f . 



(9) 



Revegf{p f , 0;e) = {g R {p } , 0; e)-gf{p f , 0; e)) tanh(e/2T) 
is the equilibrium Keldysh Green function constructed 
from the retarded and advanced ones, g^' A (p^,0;e), at 
the interface. Functions g^' (p* , 0; e) are calculated 
at the interface so the boundary condition ^ is ful- 
filled. The resulting critical current of the junctions is 
characterized by different transparency and spin-mixing 
angle and show a rich variety as seen in Figures ||. This 
quantity is defined as the maximum amplitude of the 
current reached between a phase difference of and tt. 
The sign of the critical current is cither positive giv- 
ing a "0" -junction or negative signaling a "7r" -junction. 
For arbitrary transmission and spin mixing the ana- 
lytic expression for the current is not very tractable and 
numerical analysis of the current-phase relation is more 
practical. In the two extreme limits of tunneling, t « 1, 
and high transparency, t ~ 1, the analytical expression 
is simpler and reveals the physics going on. 

Starting with the tunneling limit, the transmitted co- 
herence functions T, f as given by equations (||) are ex- 
panded to leading order in transparency T = \t\ 2 . The 
expanded T, V are then put into the expression for g R , 
eq. (^|), and to first order in T the resulting current is 
given as follows 



jtfr®)=Tj $mct> f 
J — ( 



de 

AtT 



/c + (e,e) + /c_(£,e) 



(10) 



where £±(e,0) = [ftcos(0/2) ± e sin(0/2)]- 2 tanh(^) 
and j = 2evfNfA 2 . In general the current is totally 
governed by the bound states at e = ±Acos(0/2) and 
their population at the given temperature T. Notable is 
that for all values of 0, the current-phase relation is sinu- 
soidal. Setting —JX reproduces the usual Ambegaokar- 
Baratoff expressionEJ. At = tt, is a zero-energy 
bound state and give rise to a "7r" -junction with a criti- 
cal current which increases as T _1 with decreasing tem- 
perature as shown in the inset in the left figure |[ A 
similar anomaly ria the critical current occurs for d-wave 
superconductorsEil. The difference between the anomaly 
in the two types of superconductors is that for d-wave su- 
perconductors any concentration of bulk impurities will 
give a finite width of the zero energy bound states and re- 
duce the T _1 -anomaly. For an s-wave superconductor at 
a point contact only inelastic scattering processes, phase- 
breaking impurities or, as shown below, a finite T can 
give a similar broadening and the anomaly is quite robust 
due to long inelastic (phase-breaking) scattering times, 
i.e. l/2T ine i t (phase) << A. Going away from = tt, 



£b moves to finite energy and the functions JC±(e, 0) ac- 
quire double pole structure at £{,. The two poles are 
slightly shifted and have slightly asymmetric magnitude 
in residues. Given that the residues also are of different 
sign, they contribute oppositely to the critical current. 
The separation in energy of the poles is dependent on the 
imaginary part, S, of the energy e + id. 6 can loosely be 
interpreted as an inelastic scattering rate. At T = and 
6 — > it turns out that the ground state of the junction is 
always a "0" -junction except at = tt. As temperature 
and/or the 5 is increased the region in showing a "tt v - 
state junction increases. This is clearly demonstrated in 
figure ^ where the critical current is plotted as function 
of temperature. For all but the two largest values of 
the low-T critical current is positive. 

Sticking to the tunneling limit, the results derived here 
can also be obtained using equations (93-94) of MRS. 
What is crucial to note is that in order to find the bound 
state contribution to the current the Green's functions 
must be the ones arrived at in eq. (j^). This is tak- 
ing into account the spatial dependence of the Green's 
function at the magnetic pinhole, i.e. solving the impen- 
etrable wall problem with the i'-matrix describing the 
pinhole. If this spatial dependence is neglected, and the 
bulk Green's functions are used in equations (93-94) of 
MRS, the resulting Josephson current will simply have 
two contributions sin(0 ± 0), one contribution for each 
spin band. 

Moving away from the tunneling limit the T _1 - 
anomaly is cut off by the finite value of T. Instead, the 
switching between "0" and "7r" -states happens in abrupt 
jumps in the critical current. This abruptness is spurious 
since the current-phase relation has three zeros between 
and tt and the change of state from " 7r" to " 0" occurs 
without the current being zero for every phase difference. 
Instead the junction has two local minima in energy, at 
phase difference = and tt. Right at the switching 
point the two states are degenerate and the junction state 
can be tuned with temperature. This is shown in the 
inset of the middle figure || which depicts the junction 
energy vs. phase difference for a junction with T = 0.1 
and = 37r/4. As seen, between the two energy minima 
there is a potential barrier. This barrier is at highest 
at an intermediate phase < (j> < tt where the current 
through the junction is zero. As temperature is swept 
over the switching temperature, which for this junction is 
at T sw w 0.12T C , the energy minimum jumps from <f> = tt 
to (p — as temperature is increased through T sw and 
vice versa as the temperature is decreased through T sw . 
The position in temperature depends on the two junc- 
tion parameters T and 0. For larger transparencies the 
values of tunable junctions is restricted to a decreasing 
range in just below = tt/2. 

If the transparency is taken to unity the tunability of 
the junction state with temperature vanishes. This is 
seen from the current-phase relation 
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with 

J±(0;e,e) = 



de 

4-7T 



[ft 2 



J + (<t>;e,Q)+J-(4>;e,Q) (11) 

ship ± 9) 
s 2 + A 2 cos(0±9)] tann 4T j - 



The current is now controlled by interface states located 
at £j = ±Acos[(0± 0)/2], i.e. at a position given by 
the phase 4> but shifted by ±9 for the two spin bands 



as compared to the spin-neutral case. This shows up in 
the fact that at 9 = the usual Kulik-9mel'yanchuk 
(KO) formula is recovered^ and at finite 9, eq. ( [Ti"| ) is a 
sum of two K9 supercurrents evaluated at phases shifted 
by ±9. For spin-mixing angles 9 < tt/2 the junction 
is in the "0" -state and at 9 > tt/2 in the "tt" -state. 
At 9 = 7r/2 the junction state is degenerate for every 
temperature as the current-phase relation has doubled 
periodicity. 




FIG. 3. Critical currents for different transparencies 0.01,0.1 and 0.99, left to right, and for a dense sampling of the 
spin-mixing angle, O running from to n in steps of n/20. Thick lines are in intervals of 7r/4 as a guide for the eye. All 
currents are scaled with the value of the critical current at T — and = for the current value of T. As is seen for all 
values of transparency, T, the junction may either be a "0" or a "tt" -junction depending on the degree of spin mixing. 
In the low transparency limit and at = tt the zero-energy bound state gives rise to a critical current ~ T _1 as seen in 
the inset of the left most panel. At intermediate the junction may switch between the "0" and the "tt" -junction state 
with temperature. At larger T this switching becomes more abrupt in temperature. Finally, in the high transparency 
limit the switching between "0" and the "tt" -junction state is lost and defines the junction state for all temperatures. 



V. DISCUSSION 



In this paper a general solution is derived for the equi- 
librium part of the Zaitsev-Millis-Rainer-Sauls boundary 
condition describing spin-active interfaces. This solution 
is the main result of the paper and will be an impor- 
tant part in further studies of hybrid superconductor- 
(ferro)magnetic systems. As an application, the effects 
of a magnetically active barrier, as described by a simple 
two-parameter S'-matrix, are studied. In particular, it 
is shown that spin-mixing brings about Andreev bound 
states within the superconducting gap A. The energy 
of these states is sensitive to the amount of spin-mixing 
imposed by the scattering off the interface. Comparing 
with the D9S calculated here and those obtained in the 
tunneling experiments of Stageberg et afl it is plausible 
to conclude that the spin-mixing angle is not that large 
but rather in the range |9| < tt/4 for the materials stud- 



ied in the experimenti. None the less, it is important 
to note that the shift seen in the tunneling conductance 
may be dependen±-.on, and described by, the tunneling 
barrier propertied 2 !! and thus not directly dependent of 
the exchange field in the ferromagnet probed&Ej. On ac- 
count of the Josephson coupling through a magnetically 
active interface, a small value of 9 would imply that "tt"- 
junctions are hard to realize at least in large junctions. 
To obtain a "7r" -junction it is shown that 9 must ex- 
ceed 7r/2 for any range of transparency. On the other 
hand, new experiments are in the making like.-the mag- 
netic Cobalt grains studied by S. Gueron et aO. These 
small magnetic systems may well prove to offer magnetic 
scattering where a large 9 ~ tt is realized. As an exam- 
ple, using STM-techniques, as those performed on the Au 
point contacts in Ref . [ [54| , on small magnetic grains and 
with superconducting electrodes, the Josephson physics 
described in this paper could be probed. 
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APPENDIX A: QUASICLASSICAL THEORY 

Calculations presented in this paper are done within 
the quasiclassical approximation which is a general- 
ization of the Landau Fermi-liquid theory to include 
supercondunctingta and superfluidtj phenomena. Quasi- 
classical theory is an expansion in quantities like T/Tf or 
l/£kf, which are usually of order ~ 1CP 2 — 10~ 3 in con- 
ventional superconductors. I use the quasiclassical the- 
ory for a p-rwave superfluid 3 He as worked out by Serene 
and Rainert3 together with the reaktnetal-oriented weak- 
coupling theory of Alexander et atll to describe the su- 
perconductor in proximity to a magnetically active ma- 
terial. In this appendix a brief review, or collection, of 
the building blocks of quasiclassical theory is given. 

Our starting point is the Eilenberger equation 



i VfV R g + [ef 3 



E,S] = 



(Al) 



for the 4x4 matrix propagator g(p*, R; e). Here Vf is the 
Fermi velocity, py is a point on the Fermi surface, The 
explicit 2 x 2-matrix structure of g reflects particle-hole 
(Nambu) space. The spin degree of freedom is in the pa- 
rameterization into spin scalars, g(pf, R;e), g{pf,R;e), 

f(p f ,R;e), f(p f ,R;e), and spin vectors, g(p f ,R;e), 

g(p f ,R;e), f(p f ,R;e), f(p f ,R;e) as 



9 = 



9 + 9? (f + f-<r)iv2 
«o- 2 (/-/-<t) a 2 {g - g-<r)a 2 



(A2) 



In addition to (|Al[) the propagator obeys the normaliza- 
tion condition g 2 (Pf, R', e) = — 7r 2 . There is some redun- 
dancy in the parameterization of (A2) which gives the 
following symmetries^ 



x(p f ,R;e) 



-pf,R;-s*)* 



(A3) 



where x (x) is one of the spin components g a p or / Q| g 
(ga/3 or fa/3) of the Green's function. Matsubara propa- 
gators are obtained by (e — > ie n = wrT(2n+l)), retarded 
propagators by [e — ► e + iS), and advanced propagators 
by (e — ► e — iS). Analogous symmetry relations hold for 
the self energies. 



The self energy E in equation (Al) contains impurity 
contributions, the Fermi-liquid mean fields and the order 
parameter 



E( Pf ,R;e) = Ei(p f ,R;e) + E m (p f ,R) + A(p f ,R) (A4) 

The self-consistency equations for the impurity self en- 
ergy E; given one impurity potential u\ and one impurity 
concentration n\ is 

Si (p f , R, e) = miy (p f ,p f ,R,e) 

with the quasiclassical T-matrix equation 

ii(pf,Pf,R;e) =<k(pf,Pf) 

+N f {u i (p f ,pf)g(pf, R; e)ii(p' f ' , pf , R, e))p> 

Nf is the averaged normal state density of states at the 



(A5) 



(A6) 



Fermi surface and 



denotes a Fermi-surface aver- 



age. If there are more than one type of impurity po- 
tentials equation QA5| ) will be a sum over the different 
impurity contributions, each with its own density and its 
own T-matrix. It is important to bare in mind that Ej 
is in general not diagonal in particle-hole space. The self 
energy E m contains the Fermi-liquid mean-field self en- 
ergies. It is diagonal in particle-hole space and divided 
into a symmetric (E m ) and an antisymmetric (S m ) part 
as 

Z m (p f ,R) = Tj2(A s (p f ,P f )g(p f ,R;e n )) p , (A7) 

V rn {p } ,R)=TY J {A a {pf,Pf)g{p'f,R\e n )) plf 

The Fermi-liquid interactions A^ s ' a ^ (p* , pi.) are parame- 
terized by the Fermi-liquid parameters A s and A a which 
are phenomenological parameters determined from ex- 
periments. The order parameter A is split into a singlet 
A and a triplet A parts by the singlet and triplet pairing 
interactions V 3 (pt,pt) and V t {pf,pf), and is calculated 
as 



A(p fl R) =Tj2{V s (p f ,p / f)f(p' f ,R, e n )) p , 
A(p f ,R) = TJ2 (V'iPf^fiPf, R; e n )) p , . 



(A8) 



The set of equations written above, the Eilenberger 
equation for g and the equations for the self energies 
S, must be solved self-consistently by iteration together 
with the appropriate boundary conditions imposed on 
the propagator. With g(pj,R;e) determined, physical 
quantities like the current densitymay be computed 



j(R) = 2eN f Tj2 {v f (p f )g(p f ,R;e n )) p 



(A9) 



The local density of states resolved for a given pi and a 
given spin direction e is calculated at real energies as 

N e (p f) R;e R ) = -^Im 

1 7T 



g(p f ,R;e R ) + e-g(p f ,R;e R ) 
(A10) 
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